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Abstract 



We consider a new approach to the description of dense nuclear matter 
in the framework of chirally symmetric, quark-based hadron models. As pre- 
viously in the Skyrme model, the dense environment is described in terms of 
hyperspherical cells of unit baryon number. The intrinsic curvature of these 
cells generates a new gauge interaction for the quark fields which mediates 
interactions with the ambient matter. We apply this approach to the Nambu- 
Jona-Lasinio (NJL) model, construct its curved-space quark propagator and 
solve the ladder Bcthc-Salpctcr equation for the pion. Wc find a high-density 
phase transition to chiral restoration, discuss the density dependence of the 
chiral order parameter and of the pion properties, and compare with results 
of the conventional chemical-potential approach. The new approach can ad- 
ditionally describe baryon-density-free cavities in the dense medium. 
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Typeset using REVT^ 



I. INTRODUCTION 



Nuclear matter has been the subject of intense study over the last decades It repre- 
sents in many ways the simplest nuclear many-body system and also plays significant roles 
in other areas of physics, ranging from relativistic heavy-ion collisions to the structure of 
dense stellar objects and to the evolution of the universe. In recent years, nuclear matter 
under extreme conditions, i. e. at densities a couple of times beyond the saturation density 
and at finite temperature, became a particularly active research area. 

Under these extreme conditions, the strong interactions and the hadron spectrum are 
expected to undergo qualitative changes and, in particular, transitions to new phases and 
vacua 0. Lattice simulations and different models predict, for example, chiral 
symmetry restoration and deconfinement transitions. The new phases will probably still 
be complex and nonperturbative, with an excitation spectrum containing colorless bound 
states carrying hadron quantum numbers and other collective modes, e.g. plasmons 
and plasminos Pion [§] and kaon ^ condensation and a vector-symmetric vacuum 
||10| , p!I[| have also been contemplated. Produced in finite volumes, some of the new vacua 
might furthermore have interesting coherent decay modes [|12| , p!3| . The study of all these 
phenomena promises new and direct insight into the elusive nonperturbative sector of the 
underlying theory, QCD. 

Since the advent of relativistic heavy-ion accelerators at Brookhaven and CERN in 1986, 
matter under extreme conditions became also accessible in the laboratory. In central heavy- 
ion collisions very dense nuclear matter, albeit probably not close to equilibrium, can be 
produced. Fixed-target experiments with heavy beams, as e.g. the now operational gold 
beam at the Brookhaven AGS and the planned lead beam at CERN's SPS, will produce 
particularly large final state multiplicities of a few thousand charged particles per central 
collision, and accordingly very large baryo-chemical potentials. 

Realistic calculations of nuclear matter properties at these high densities are difficult 
to perform. Direct information from the lattice cannot be expected soon, since a chemical 
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potential renders the euclidean action of QCD complex and conventional Monte-Carlo tech- 
niques inapplicable [|I^. The traditional approaches, as for example Brueckner-Hartree- 
Fock and variational calculations, are based on effective interactions {e.g. from meson- 



exchange models [^) between point-like nucleons. At high densities, however, this neglect 
of the nucleon structure is a serious shortcoming, since the chemical potential in compressed 
nuclear matter is of the order of or larger than the lowest nucleon excitation energies. 

Additionally, density-induced modifications of the hadron structure, which recently be- 
came an active field of study |jl6|, should be taken into account and will in turn affect the 



matter properties. Particularly dramatic changes in the hadron structure can be expected 
in the vicinity of phase transitions, where the ground state of the matter rearranges itself in 
a fundamental way. The chiral phase transition alters e.g. the hadron spectrum drastically, 
and the deconfinement transition transforms hadrons either into weakly bound states or 
dissolves them completely into a quark-gluon plasma. 

To account for the hadron structure and its modifications in high-density nuclear matter 
calculations is a challenging task. Indeed, a realistic treatment is beyond the reach of present 
capabilities. One can, however, attempt a more qualitative description in the framework of 
baryon models, based on the unit-cell method. This type of mean-field approach is familiar 
in solid state physics for the description of homogeneous or periodic media and becomes 
reliable at very high densities. Also, it is simple enough to be applicable in many existing 
baryon models. 

Let us briefly recapitulate the basic concepts underlying the unit cell approach, which 
we will later generalize. After dividing the nuclear matter into identical cells of unit baryon 
number, periodic (or "twisted" periodic) boundary conditions are imposed on the cell, i.e. 
the physics on opposite surfaces is identified. As a consequence the cell looses its boundary 
and becomes self-contained, while the presence of the neighbors is reflected in its periodic 
images. In topological terms, the cell becomes a three-torus x x . 

The size and shape of the unit cell is found by minimizing its energy under variations of 
its geometry. In this way interactions with its surroundings determine via the variational 



principle the detailed cell structure. This is exactly how one derives, e.g., the simple cubic 
unit cell of sodium chloride from the Coulomb interaction. Information about the neighbors 
and the mutual interactions is then encoded in one single, isolated cell. In the context of some 
baryon models, and soliton models in particular, the unit cell approach has the additional 
advantage of describing the density- dependent structure of the baryons and their mutual 
interactions in matter consistently, namely on the basis of the same underlying dynamics. 
In the Skyrme model - a widely studied soliton model based on an effective chiral 



meson lagrangian - a face-centered cubic (fee) unit cell was found along these lines a few 
years ago [|l^ . More recently, Manton and Ruback pointed out that the cell energy could be 
lowered further by giving the unit cell an intrinsic curvature ||T9|j20[| . Their proposal can be 



regarded as an extension of the variational principle, which now allows all the geometrical 
properties, not only the shape, of the cell to be determined by the given dynamics. The 
optimal cell turned out to be a three-dimensional hypersphere S^{L) of radius L. Exclusively 
in this geometry the skyrmion can attain its absolute energy minimum, the Bogomolnyi 
bound |T§. 

The simple form of the optimal unit cell, S^, is essentially determined by chiral sym- 
metry, is invariant under the group SO {4), which is (locally) equivalent to the chiral 
group SU{2) X SU{2). The pion fields therefore take values on S"^, and their lowest-energy 
configuration is attained if they are defined on another S*^, since the resulting map minimizes 
the deformation energy stored in the field gradients. 

The most remarkable new feature of hyperspherical cells is their curvature, which induces 
additional dynamics for the fields. According to the unit cell concept, these interactions are 
ascribed to the dense environment. As in conventional unit cells, the cell boundary has 
disappeared. The periodicity in the three angular coordinates of can be regarded as a 
reflection of the neighboring, identical cells. The replacement of the surrounding cells and 
their boundaries by periodic coordinates is, as pointed out before, a generic feature of most 
unit cells and it supports the physical intuition in dealing with the hypersphere. 

The hypersphere approach has by now been extensively applied in the Skyrme model. 



Perhaps most remarkably, a chiral phase transition emerged at high densities, accompanied 
by parity doubhng of the hadron spectrum and the disappearance of the Goldstone modes 
||2T|. Also, strangeness condensates with distinct characteristics due to changes in the ex- 



tended baryon structure have been found |22], and the vector limit ||T0| phase transition 



has been studied |TT||. Furthermore, bag formation and dissolution some properties of 



the loop expansion in the presence of a skyrmion p4| and its excitation spectrum |25| have 
been investigated. Skyrmion-like profiles on the hypersphere have also been derived from 
instantons |2^, following an idea by Atiyah and Manton ||2^, and other soliton solutions 



have been constructed 28 



A remarkable and generic finding of these studies is that hypersphere calculations repro- 
duce almost quantitatively the results of similar fee array calculations, whenever those are 
available. This holds in particular for the chiral phase transition, its critical density, the den- 
sity dependence of the order parameter and the energy density, and it established confidence 
in the new approach. At the same time the calculational effort is greatly reduced. Array 
calculations require the numerical solution of partial-differential field equations on a three- 
dimensional grid. Hypersphere calculations, on the other hand, involve at most the solution 
of an ordinary differential equation, as for isolated skyrmions. Some results can even be 
obtained analytically and allow a straightforward and transparent physical interpretation. 

The rather successful and consistent results of this unorthodox approach call for both a 
better understanding of its physical foundations and for further exploration of its range of 
applicability. In the present paper we will take a first step in this direction by extending 
the hypersphere approach to quark-based hadron models, and to the Nambu-Jona-Lasinio 
(NJL) model ||2^ in particular. Some results in the chiral symmetry breaking sector have 



already been published in ref. pO| . 

The NJL model shares with the Skyrme model the underlying chiral symmetry, which 
led to the geometry of 5''^. In many other respects, however, the NJL dynamics is different. 
In contrast to the purely mesonic degrees of freedom in Skyrme-like models, it is based solely 
on quark fields. Furthermore, it does not describe nucleons as solitons and it breaks chiral 



symmetry dynamically. These differences can help to disentangle effects specific to the 
Skyrme model from more generic or even model-independent features of the hypersphere 
approach. In particular, they could help to clarify whether the specific geometrical and 
topological properties of the Skyrme model are indispensable for the effectiveness of the 
hypersphere approach. 

The choice of unit cells even beyond the Skyrme model is additionally motivated 
by their maximal symmetry under all three-dimensional curved spaces. Hyperspherical 
cells can thus be considered as the simplest generalization of flat unit cells, and many of 
their characteristic features are indeed reminiscent of flat cells. This tendency to preserve 
quahtative properties of euclidean space is enhanced by the conformal equivalence between 

and R^. The use of cells furthermore avoids potential problems with explicit chiral 
symmetry breaking by a cell boundary (as, for example, in bag models). 

An alternative, more conventional approach to dense matter in the NJL model, in which 
the quarks are coupled to a homogeneous, external baryon number source, has been studied 
by Bernard et al. [Q and by others [^. As in the Skyrme model, where the comparison 



with conventional array calculations served as a test for the hypersphere approach, we will 
proflt from comparing our results to those of the chemical potential approach. 

Let us add a couple of comments concerning the uses and limitations of our approach. 
Any attempt to describe nuclear matter in terms of unit cells, i. e. in a mean fleld framework, 
is bound to be qualitative at best . Approximating the short- and medium-range nucleon- 
nucleon correlations by a strict long-range order and neglecting the nucleons kinetic energy 
are of course rather severe simpliflcations, except at very high densities such as in neutron 
stars, where nuclear matter might even crystallize. 

On the other hand, one may expect that bulk features of nuclear matter are robust enough 
to survive these approximations, and that they can be revealed in the unit-cell framework. 
Global features and driving mechanisms of high-density phase transitions, for example, could 
fall into this category, and the alternative Fermi-gas treatment will fail to describe them. 
As already mentioned, the unit cell approach offers the so far unique possibility to account 



consistently for the extended baryon structure, which certainly contributes qualitatively 
important physics at high densities. 

The paper is organized as follows: In section II we establish the general framework for 
our investigation by implementing the quark fields into the curved unit cell. Section III deals 
with the specific NJL dynamics and the construction of its curved-cell quark propagator in 
Hartree-Fock approximation. The latter is a key step in our program. Its more technical 
parts are relegated to appendices A and B. 

In section IV the chiral symmetry breaking sector and its density dependence is ex- 
amined. We derive, in particular, the gap equation, which describes the dynamical quark 
mass generation, and the expression for the chiral order parameter. In section V the Bethe- 
Salpeter equation in the quark- antiquark channel is set up in the curved background, and 
its solution for the pion in ladder approximation is found. 

Prom the explicit pion wave function we calculate in section VI the pion decay constant 
and its density dependence. In section VII the zero-density limits of the various calculated 
quantities are derived and compared with the standard flat-space NJL results. Section 
VIII discusses the quantitative aspects of our results and compares them with results of 
the chemical potential approach. Finally, section IX summarizes the paper and offers some 
conclusions and perspectives for future research. 

II. QUARKS IN THE CURVED UNIT CELL 

In this section we describe the implementation of quark fields into a curved unit cell and 
summarize the generic physical implications. This discussion is quite general and provides 
a framework for the adaptation of a large class of quark-based baryon models. 

As motivated in the introduction, our unit cell has the geometry of a three-dimensional 
hypersphere, S^{L), of radius L. Accordingly, the flat Minkowski metric is replaced by the 
metric of S^{L) x R {R indicates the real, flat time dimension), 

ds^ = g^.^dx^'dx'' = de - L^[dn^ + sin^ i_i{de^ + sin^ edc^^)], (1) 

8 



which we write in the polar coordinates of R^, i.e. {x'^} = {t, /i, 6, 0). Its intrinsic curvature 
leads to new interactions for the spinor fields, as it does, for example, in general relativity. 
In our context, however, these interactions will be ascribed to the dense environment. 

The spatial section of (0) has six symmetry generators^, the maximal number in a three- 
dimensional, curved space. Three of them are just the ordinary rotations in R^, while the 
remaining three are compact analogs of Lorentz boost generators and can be regarded as 
"generalized translations" in 5^. Combined with the usual translations in the flat time direc- 
tion, the metric (|l|) therefore possesses seven symmetries, only three less than the Minkowski 
metric of R^. This large symmetry group will be helpful both in establishing analogies with 
fiat unit cells and in explicit calculations. Together with the conformal relation between the 
hypersphere and Minkowski space, it often leads to structural similarities in expressions for 
physical quantities. 

In the Skyrme model, which contains only spin-0 pion fields, the change of the lagrangian 
to the new metric was sufficient to adapt the model to the hypersphere. Spinor fields in 
curved space, however, require the introduction of additional concepts. The straightforward 
generalization of the fields from representations of the Lorentz group to representations of 
the general linear group of (infinitesimal) coordinate transformations fails for fermions, since 
the linear group has no spinor representations. 



This impasse can be circumvented by defining the spinor fields in a local orthonormal 



basis, given by the vierbein fields e"(x). The e"" are generated by linear, space-time dependent 
transformations of the coordinate basis, written in terms of the 4x4 matrix e°^(x): 

e^ix) = e''^{x)dx'', a €{0,1,2,3}. (2) 

The corresponding vector fields, which form the components of the gradient operator in the 
orthonormal basis, can be obtained from the inverse vierbein coefficients e^'^, defined by 



^ These SO (A) generators form together with the time translation generator the Killing vectors of 
the metric (||). 

9 



e\{x)e,^{x) = 5t. (3) 

The Cj,'^ transform coordinate-vector components into Lorentz-vector components (in the 
local frame) and are obtained from the e"^^ by lowering the Latin index with the Minkowski 
metric {rjah] = diag{l, —1, —1, —1) of the local frames and by raising the Greek index with 
the hypersphere metric g. The gradients are then given by 

eaix) = e/(x) d,. (4) 

Due to its orthonormality, the vierbein is simply related to the metric (|l]): 

9^.u{x) = Vab e"^(x) e\{x). (5) 

In our context it is crucial to note that this decomposition of a given metric in terms of 
the vierbein frames is not unique. Indeed, an arbitrary, local Lorentz transformation of 
the Latin vierbein indices leaves the expression unchanged and thus leads to the same 
metric. These Lorentz rotations consequently form a gauge group, which just changes the 
relative orientation convention (the "gauge") of the local frames, but leaves the physics (i.e. 
the metric) invariant. It generalizes the global Lorentz group of Minkowski space, and its 
spinor representations allow the definition of fermion fields in curved space. In order to 
preserve the local Lorentz invariance of the lagrangian, gradients of the quark fields have to 
be replaced by gauge-covariant derivatives, 

ej'{d^ + ^i^)(l= {ea + ^a)q, (6) 

which transform homogeneously under the gauge group. In eq. (^) we have introduced a 
(4 X 4-matrix-valued) gauge field, the spin connection 

uJa = e/ uj^ = --uj^'f^bc (7) 

{o'ab = f [7a,7f)]), which just compensates the change of the spinor field components due to 
different frame orientations at neighboring points by a spatial 5*0 (3) rotation. In this paper, 
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we take uo to be the Levi-Civita connection^, which is uniquely determined by the metric 



and the vierbein p3| : 

(5^9^-r/'')e/. (8) 



^ fee ^ _g6 gC 



The Christoffel connection F is the corresponding gauge field of the local Lorentz group in 
the vector representation and can be directly obtained from the metric: 

r/" = \g.,{d''g'"^ + dPg'^'' - d'^g^"). (9) 

The explicit form of uj depends, of course, on the gauge. The group structure of S"^ 
implies the existence of a particularly useful choice, which we will refer to as the "Maurer- 
Cartan gauge". In it, the spin connection takes a simple and, in particular, space-time 
independent form: 

a;a = — eafeccr*'' (a, 6, c e {1, 2, 3}) and = 0. (10) 

[pah = f [7a, 7fe]) More details about this gauge, in which all of the following calculations will 
be performed, can be found in appendix A. 

In the context of nuclear matter unit-cells, it might at first be tempting to think of the 
spin connection as a covariantized chemical potential, since it acts as a constant, isoscalar 
vector interaction for the quarks and grows with density. Such an interpretation would, how- 
ever, be misleading. Besides its gauge dependence, oj has a different Dirac-matrix structure 
and its time component vanishes identically in the matter rest frame. 

Following the procedure outlined above, one can now adopt a general fermionic lagrangian 
from Minkowski space to the hypersphere. The corresponding action is obtained by inte- 
grating the lagrangian over S^iV) x i?, i.e. 



^The Levi-Civita connection is metric and torsion-free (see appendix A) and thus corresponds to 
the natural and minimal generalization of flat unit cells. In principle, one could use a more general 
spin connection, which would introduce a non-zero torsion into the unit cell. This issue might 
deserve future investigation. 
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S = j d^{x) C{x). (11) 

The integration measure dfi{x) contains the jacobian a/— g, where g(x) = — sin^ /x sin^ ^ 
is the determinant of the metric tensor (|^): 

djj,{x) = \J —g{x) d^x = sin^ sin 6 djj, dO dcf) dt. (12) 

In order to complete the definition of the unit cell it remains to set its total baryon 
number to one. Since this is most conveniently done at the level of the quark propagator, 
we will reserve this last step for the next section. 

Let us now summarize the model-independent physical implications of the curved unit 
cell for fermions. The novel aspects can be classified according to their either local or global 
nature. The local effects are due to the finite curvature: gradients are replaced by vierbein 
vector fields and a new gauge interaction emerges. 

However, as already pointed out in the introduction, there are also new features of global, 
i.e. topological character. Due to the compactness of the hypersphere, the fermion spectrum 
(see appendix A for its explicit form) will be discrete, with the energy scale of the excitations 
determined by the inverse radius of the cell. Furthermore, potentially symmetry-breaking 
cell boundaries are absent. The periodicity in the three (angular) spatial directions fi, 6 and 
(j) is reminiscent of the periodic boundary conditions in fiat unit cells and can be regarded 
as a refiection of the neighboring, identical cells. Quarks moving around the cell more than 
once are accordingly interpreted as entering a neighboring cell. This physical interpretation 
will become explicit in particular during the construction of the fermion propagator. 

We conclude this section with a remark on the perturbative treatment of space curva- 
ture, which is frequently used in general relativity. While this approach often significantly 
simplifies exphcit calculations, it will be insufficient for our purpose, for two reasons: First, 
the curvature of our unit cell is of the order of V^~^, where V is the average volume occupied 
by one nucleon in nuclear matter, and thus becomes large at high densities. Secondly, a 
perturbative approximation to the metric disregards the global structure of the unit cell. 
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This is another serious shortcoming at high densities, where the cell becomes smaller and, as 
we will show in the next section, the quarks get lighter on the approach to chiral restoration. 
Under these circumstances the global features of the cell become important. 



III. NJL MODEL AND QUARK PROPAGATOR ON S'-^ x R 

Up to now our discussion established a generic framework for the study of fermionic 
models in hyperspherical unit cells. In order to proceed further and to arrive at quantitative 
results, we have to specify the dynamics of the quarks. 

Many of the existing, chirally symmetric quark and quark-meson models are accessible 
in our approach. The Nambu-Jona-Lasinio (NJL) model ||29|, however, seems particularly 



well suited for a first, explorative study aimed at testing the use of hyperspherical unit cells 
beyond the Skyrme model. As already indicated in the introduction, the NJL dynamics is 
simple and contrasts the Skyrme model in many aspects, except for their common chiral 
symmetry. For example, the NJL model is entirely quark-based (in its modern form), it is 
not a soliton modelQ and it breaks chiral symmetry dynamically. 



In fiat Minkowski space, the fiavor-SU (2) version of the NJL model is based on the 
lagrangian 

^=2 [^^^^i^^ - i9^,q)Y(l] - rrioqq + g (gL^g) (gL^g). (13) 

The g's denote Dirac-spinor quark fields, which are color triplets and isospin doublets. The 
scalar and pseudoscalar four-quark contact interactions resemble instanton-generated ver- 
tices of light quarks in QCD and are written in a shorthand notation with the help of 



^In this paper we will consider only the original, fermionic version of the NJL model. A restricted 
bosonization of some generalized NJL lagrangians can support stable, skyrmion-like soliton solu- 
tions |34|. Their eventual investigation on could shed more light on connections between the 
hypersphere results in the Skyrme and NJL models. 
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the four matrices r° = l^lr, = 275X4 {i = 1, 2, 3). (The index a G {0, 1, 2, 3} is imphcitly 
summed over in the interaction term above.) In most of the following discussion we will set 
the small current quark masses rrio to zero. 

Due to its dynamical symmetry breaking mechanism, simplicity and phenomenological 
successes, the NJL model is widely employed as an effective low-energy theory for QCD 
36| , |37|| . Recently, the original version of the NJL model has been generalized in a couple of 



directions, e.g. by adding vector and axial-vector couplings or by extensions to the SU{3)- 
fiavor sector. These developments and many applications to meson and baryon physics, 
including some in hot and dense nuclear matter, are reviewed in refs. [^,0. We will 
work with the original, "minimal" version of the model, however, since our aim is not a 
detailed phenomenological analysis, but rather the investigation of generic aspects of the 
new approach. 

We begin our study by adapting the NJL lagrangian ([T3|) according to the procedure 
established in the preceding section to the hyperspherical unit cell, i.e. we replace the flat 
Minkowski metric by eq. (|I|) and the gradients by gauge- covariant derivatives: 

^ = ^[Ql''{ea + uJa)q-q{'ea-uJah''q] - moqq + g (qVaq) (qV^q). (14) 

(The arrow indicates that the corresponding differential operator acts to the left.) 

Most of our following analysis will be based on the standard Hartree-Fock approximation 
to the NJL dynamics p9[. As is well known, the quarks dynamically acquire already at 



this level a finite self-energy, if the coupling g exceeds a critical value. Due to the pointlike 
character of the interaction, the self-energy is space-time independent. The additional gauge 
interaction on the hypersphere preserves this property, since the spin connection is a non- 
propagating background field. In the presence of a finite baryon density, the self-energy 
acquires besides the standard scalar piece, which acts as a "constituent" mass, an additional 
vector part: 

E = - 7oS„. (15) 
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All the vacuum, quark and pion properties of interest in this paper can be conveniently 
calculated with the help of the NJL constituent quark propagator, which includes the self- 
energy (|15|). The next and pivotal step of our program is therefore the derivation of an 
explicit expression for this propagator on x R. 

To this end, and at the densities of interest the crucial modifications in the quark dy- 
namics due to the large curvature have to be taken into account exactly. While in general 
analytical expressions for propagators in curved spaces can rarely be found, the high sym- 
metry of the hypersphere allows us to derive the quark propagator on x R exactly and 
in closed form. Our derivation starts from the definition of the propagator as the inverse of 
the gauge-covariant Dirac operator, 

[«7"(ea + ^a) + Klo - {mo + S,)] S{x,y) = S^{x,y). (16) 

(The explicit form of the delta function on 5*^ x i? is given in appendix B.) We then specialize 
to the Maurer-Cartan gauge, separate S into invariant scalar amplitudes and rewrite them 
such that the dependence on all but the geodesic coordinate of is removed. The remaining 
/i dependence can then be absorbed into redefined amplitudes, which finally allows the 
inversion to be performed by Fourier methods. The detailed calculation is described in 
appendix B. 

The resulting Hartree-Fock propagator on x R, in the chiral 7-matrix basis, reads 

Six) = Six, 0) = e^(^^^+^"*) [ Soifi, t) 70 - S.ifi, t) n - S2ifi, t)]. (17) 

The exponential factor, which contains the Dirac spin matrix^ S, is the spin parallel prop- 
agator and fi6, 0) denotes the usual unit-vector of euclidean R^ in spherical coordinates. 
The decomposition of the propagator into the three amplitudes 



is the direct product of the 2x2 unit matrix and the standard Pauh matrices, see appendix 

B. 
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So{^i,t) = -^a{sm^^)-' ^ (-1)"[2/; + tan ^J„], (18) 

n=—oo 
-{-oo 

5i(/x,t)=z^(sin/.)-i E Mr [2/:- cot ^/;], (19) 

^ n=-oo ^ 

52(/i,t) =ma(sm/i)-i E (-ir[24 + tan^/n], (20) 

n=— oo 

(a = l/47rL^) reflects the symmetry properties of 5*^ x R. Tlie Si{fi,t) are expressed as 
sums over all geodesic patlis on wliicli a quark can propagate between tfie pole {y = 0) 
and tlie point x. Since tlie cell is compact, these paths contain n full circles around the 
hypersphere. This becomes transparent by noting that all the integrals /„ = /(/i + 2mT,t) 
and their derivatives with respect to fi (t), denoted by primes (dots), are derived from the 
two-dimensional propagator 

dko r dk Qi{ktiL-kot) 



(m = mo + Ss). The additional /i dependence in equations ([T8|) - (^Of) just translates 
propagation along a fiat "/i-axis" into geodesic propagation in the /i direction on the curved 
hypersphere. 

The propagator (|17D describes quark propagation in a cell with zero net baryon number 
and consequently satisfies the usual Feynman vacuum boundary conditions. Its adaptation 
to the cell's physical ground state with baryon number one is, however, straightforward. 



owing to the transparent implementation of the boundary conditions in eq. (pl|). 

The "Fermi sea" of three valence quarks in the cell affects the singularities of the inte- 
grand of J(/i, t) inside the Fermi sphere both by the appearance of poles from positive energy 
holes and by the suppression of poles from Pauli-blocked states. Accordingly, the integrand 
in eq. (^Tf ) has to be modified to 

^i{kfiL—kot) 



e 



1 777 

—6{ko - uj)Q{kF - \k\) 



(22) 



kg — uj'^ + ie uj 

where u = y/k"^ + m? and kp is the Fermi momentum of the quarks. With the explicit 
expression for the propagator at hand, we are now ready to enter the discussion of the quark 
dynamics in the hyperspherical cell. 
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IV. DYNAMICAL CHIRAL SYMMETRY BREAKING 



In the Skyrme model, the perhaps most important and unexpected result of the hyper- 
sphere approach was the prediction of a chiral restoration phase transition at high densities 
^,21|. However, the restoration mechanism exploits specific topological and geometrical 
features and, in particular, the hedgehog structure of the skyrmion]^. It is therefore not a 
priori clear if and how hyperspherical unit cells can give rise to chiral restoration in other 
models. 

In the NJL model, in particular, chiral symmetry breaking occurs dynamically, by quark- 



antiquark pair condensation in the vacuum p9[. At the same time the quarks become 



"dressed" by a virtual quark-antiquark cloud and acquire a "constituent" mass. This mech- 
anism, analogous to the BCS mechanism for gap formation in superconductivity and perhaps 
mediated in a similar way by instantons in QCD |3^, is fundamentally different from its 



Skyrme model counterpart. In the present section we will examine how it is affected by a 
finite baryon density, described in terms of hyperspherical unit cells. 

We base our study, following Nambu and Jona-Lasinio, on the Schwinger-Dyson equation 

E{x,y) = S,\x,y)-S-\x,y) (23) 

for the quark self-energy, since the development of a constituent quark mass is a direct 
signature for chiral symmetry breaking. Equation (|2^) has the same form as in fiat space 



and its derivation proceeds along the same standard lines 39 . Note, however, that both S 



^In the Skyrme model in flat space, chiral symmetry is broken by a nonlinear parametrization of 
the pion fields. On the hypersphere at high densities the ground state changes into a generalized 
hedgehog form. In addition to the coupling of isospin transformations and rotations of the standard 
hedgehog, it also couples the axial generators of the chiral group to the "translations" on the 
hypersphere. Projection onto the physical rotation and translation eigenstates then reveals the 
restoration of the full chiral symmetry [^] . 
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and its free {i.e. g = 0) counterpart 5*0 contain the interactions with the curved background 
and, in particular, with the gauge field u). 



It is well known that nontrivial solutions of eq. (|23D, i.e. finite constituent masses, can 
only be generated nonperturbatively. Both the Hartree approximation, which becomes exact 
in the large-rzc limit of the NJL model (n^ is the number of colors of the quark fields), and 
the Hartree-Fock approximation are frequently used for this purpose. In the following, we 
will adopt the Hartree-Fock approximation to allow for a direct comparison with results of 
the conventional approach [Q. To find the Hartree-Fock solution of eq. (^), we start from 
the solution to 0{g), i.e. to first order in the NJL interaction: 

E(x,i/) = S(x)5^(x,y), (24) 

with 

S(x) = S = ig{2Tatr [S'o"(x,x)rj - [S'^f (x,x) + Sq{x,x)] F^}. (25) 

Due to the pointlike character of the NJL interaction the self-energy is, as anticipated, 
space-time independent. For the same reason, both coincidence limits of the propagator 
(p!7|) appear above, 

S^{x,x) = 5±(0) = hm 5(xo,x = 0) = 4±7o - ^2, (26) 



and are independent of x. Note that the amplitude Si in eq. (jl^) has to vanish in the 
coincidence limit, since it multiplies the unit distance vector. 

Both 5*0 and 5*2 contain (quadratic and logarithmic) short distances singularities, which 
originate from the integral I{fi,t), eq. (|21|), and its derivatives. We therefore multiply the 
corresponding integrands by a smoothed theta function^ 

e.(A-|A;|) = ^,,(l,l_\)^^ , (27) 



^Due to its negative mass dimension, the NJL coupling is not renormalizable. The introduced 
cutoff A is therefore physical [Q. 
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which damps the contributions from the high-momentum region and thus regularizes the 
integrals. Note that cutting off the spatial part of the momentum integration is consistent 
with the symmetry properties of the metric and of the "Fermi sea" of valence quarks. 

The representation of the propagator amplitudes as a sum over paths, eqs. ( pTSD - (pO[), 
is sometimes inconvenient for practical purposes. At high densities, in particular, the unit 
cell radius L is small and a large number of paths contribute. In appendix C, these sums 
are therefore transformed into equivalent mode sums over the discrete Dirac spectrum on 

X i?, with the energies 

2n + l 



ujn = ^kl + w?, kn = ^j^, riG {1,2,3,...} (28) 

and the corresponding degeneracies 

Dn = 2n{n + 1) (29) 

of the nth quark level (for a given flavor and color), which emerges in the course of the trans- 
formation. (An alternative, more direct derivation of the spectrum is given in appendix A.) 
In the spectral representation, the regularized coincidence limits of the propagator functions 
take the form (cf. appendix C) 

i °° 

So = Tj77 E ^ne,(A - knMkp ' kn), (30) 

'^^ 71=1 

im °° D 

S2 = j^Y: , " e.(A - kr.)Q{k^ - kp), (31) 
41/ n=i Jkl +im? 



where V = 27r^L^ is the volume of the unit cell. 

From the Dirac spectrum we can immediately read off the Fermi energy up, since all 
three valence quarks are contained in the lowest energy level: 



ujp = ^kj, + m2, kp = ki = — . (32) 

Since this level is only partially flUed, we have to prevent its remaining states from being 
counted in the above mode sums. To this end we associate a flUing factor 



19 



with each level, where Nn is the number of valence quarks in the nth level, i.e. /i = | and 
/„ = for > 1. 

We are now ready to give explicit expressions for the self-energy, which properly take 
the valence quarks into account. Decomposing the Dirac- matrix structure of the self-energy 
in the rest frame into scalar and vector parts, 

S = - 7oS„, (34) 

as suggested by the symmetry properties and anticipated in the last section, and inserting 
the coincidence limit of the propagator into eq. (|25|) , we obtain 

= -4:ig (1 + 2ncnf) §2 

= ^(1 + 2n^nf) £ ^"(^-^») e,(A - (35) 

V n=l 

t DJ^ = (36) 



V V 

n=l 

The above expressions exhibit clear analogies to their flat-space counterparts and allow a 
direct physical interpretation. The scalar self-energy receives contributions from virtual 
quark loops in all accessible levels of the Dirac sea down to the smooth cutoff A. The three 
occupied states in the Fermi sea contribute with opposite sign. In both self-energies the 
finite-density effects, mediated by the curved cell, manifest themselves mainly through the 
altered quark momenta kn and degeneracies which replace the plane wave spectrum of 
fiat space. 

The vector part of the self-energy is entirely due to the three valence quarks and vanishes 
in an "empty" cell, even at finite ambient baryon density. This is in contrast to the chemical- 
potential approach of ref. [^], where is necessarily finite at finite density. Furthermore, 
S^, does not depend on the value of the Fermi energy, but only on its density-independent 
relative position in the spectrum. 
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We note in passing that the not exphcitly rzc- dependent terms in the scalar self-energy 
(p5|) arise from exchange diagrams and become negligible in the large-ric limit of the NJL 
modelQ. They are kept, however, in the Hartree-Fock approximation. 

Since a constant self-energy just shifts the mass and energy scales of the quarks, the 
0{g) expressions in eqs. ( P^D and (^) have already the functional form of the Hartree-Fock 
solutions. We still have to impose the self-consistency condition, however, which requires 
the mass in the quark loop-propagator on the right-hand side of eq. ( p5|) to be identical to 
the scalar part of the self-energy itself^ and leads to the gap equation 

m = ^(1 + 2n,nf) £ ^"(^-^-) e,(A - k^). (37) 

V ^ Un 

At least at low densities and for sufficiently large coupling g, we expect the perturbative 
vacuum to become unstable, as it does in flat space. The simultaneously generated quark 
mass can then be found as the nontrivial solution of eq. (^). The additional, trivial solution 
m = does of course always exist. 

With the self-consistent self-energies at hand, the Hartree-Fock propagator is now com- 
pletely determined and vacuum expectation values of bilinear quark operators can be cal- 
culated. In the remainder of this section, we will specifically consider the vector and scalar 
quark condensates. The vector condensate, or equivalently the quark number density, pro- 
vides a consistency check on our treatment of the valence quark sector, whereas the scalar 
condensate plays a central role in the discussion of dynamical chiral symmetry breaking and 
its density dependence. 



^At the physical nc = 3 the exchange contributions are an order of magnitude smaller than the 
direct terms. 

^Since our focus is on spontaneous chiral symmetry breaking, we neglect the small explicit break- 
ing and specialize to the chiral limit, mo = 0. The generalization to a finite current mass is 
straightforward. 
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Let us start with the vector condensates. They are equal to the quark number densities 
in the cell, which add up to 



< u^u > + < (fd >= — , 



(38) 



corresponding to three valence quarks per cell. We further specialize to isosymmetric nuclear 
matter, so that the quark number densities of both up and down quarks become equal: 



< u^u >=< d^d >=< q^q >= — . 

V 



(39) 



To check our treatment of the valence quark sector, we now calculate < q^q> for comparison 
directly from the coincidence limit of the quark propagator, eq. (p6|): 



<q'q> - ov.^^c 



— itr. 
—iti. 



7,c 



7o(5-(0)p-S(0),=o^ 



7o5oJ =-/^i/i = — . 



(40) 



As expected, the two results agree and thus confirm the correct pole structure of the prop- 
agator. 

We now turn to the scalar quark condensate, which is the standard order parameter of 
the chiral phase transition. Again, it can be directly obtained from the propagator: 



< qq > 



linitr^,c < Tg(x = 0,t)g(0) >= -itr^,^ [^"(0) 



Tlrm 



V 



n n _ f 



n=l 



(41) 



Note that < gg > is proportional to the quark mass, i.e. to the solution of the gap equation 
(^). Chiral symmetry breaking manifests itself, as anticipated, simultaneously in a finite 
quark condensate and a dynamically generated quark mass. If a nontrivial solution of the 
gap equation exists and if it is energetically favorable compared to the other solutions (see 
the discussion in section VIII), a quark condensate necessarily develops. 

The gap equation (|37| ) can be solved numerically and indeed yields, at low densities and 
for sufficiently large couplings, a finite quark mass. As expected, the mass decreases with 
density and finally vanishes, signaling the chiral restoration transition. The quantitative 
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solution of the gap equation as a function of density and the evaluation of the associated 
free energy will be subject of section VIII, where we also discuss the behavior of the scalar 
quark condensate and the order of the chiral phase transition. Section VIII furthermore 
contains a comparison of our results with those of ref. 0]. 



V. THE PION WAVE FUNCTION 

As a consequence of Goldstone's theorem, the spontaneous breakdown of chiral symmetry 
is accompanied by the appearance of an iso-triplet of (almost) massless pions. They emerge 
in the NJL model as quark-antiquark pairs, tightly bound by the same interaction which 
leads to quark condensation. In the present section we study this mechanism and its density 
dependence in our framework. 

Our strategy will be to generalize the original NJL approach of calculating the pion wave 
function to the curved unit cell. Most of the new features encountered are related to the 
changes in both the quark and pion spectra. (For the actual spectra of free fermions and 
spin-0 bosons on x R see appendices A and D.) As a consequence, the Fourier transform 
of the spatial coordinate dependence becomes impractical and our calculation will be done 
in coordinate space. 

Following NJL ||29|, we start from the connected four-point function in the quark- 
antiquark channel, 

S'^^\xi,yi;y2,X2)ai3;S'r =< qa{xi)qi3{yi)q^{x2)q5{y2)\0 >c (42) 

(flavor and color indices of the quarks are suppressed), with the aim to derive an equa- 
tion for its pion pole contribution. As in flat space, the Green function (H2) satisfies an 



inhomogeneous Bethe-Salpeter equation PT| , 



S'-^\xi,yi;y2,X2)ai3;5'r = Go{xi,yi;y2,X2)ai3;S'r + J d^{zi) j dfi{z2) J dfi^zs) J dfi{z4) 

, Z2', Z3, Z4)a'l3';5'Y ,Z4;y2,X2);5'YS'y (43) 
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(the integration measure d^{z) was defined in eq. (|T2D), which can be derived either by 
summing the rescattering series of interactions specified by the gg-irreducible vertex kernel 
or formally from the generating functional of the NJL model. In eq. (^) we introduced 
the quark- ant iquark propagator 

Go{.xi,yi,y2,X2)ap-5-^ = S{xi,y2)asS{x2,yi%^, (44) 

where S is the quark propagator on x R. 

Mesonic bound states appear as poles in eq. (^3|) . In flat space one finds the correspond- 
ing intermediate state wave functions after a four-dimensional Fourier transform as the pole 



residua, which satisfy a homogeneous Bethe-Salpeter equation In general curved space- 
times the situation becomes more involved. Due to the lack of translational invariance (and 
consequently the absence of a natural global coordinate system, as provided by the cartesian 
coordinates in fiat space), the Fourier transform would have to be replaced by a generalized 
spectral transform in accord with the symmetry properties of the metric^ 

Fortunately, however, the time direction of our metric is fiat and we can isolate the 
pole piece by a Fourier transform of the time dependence only. Accordingly, we write the 
Bethe-Salpeter amplitude of the pion as 

< 0|rg„(xi)g^(2/i)|7r°(p) >= X°(a^i, j/i; p)a/3, (45) 

J2ujpV 

where p = {ujp, p} stands for the set of four quantum numbers which replace the fiat-space 
four-momentum vector of a free pion. The pion energy Up and the conserved "momentum" 
quantum numbers, p = {n,l,m}, specify the spatial state of the pion completely and are 
discussed in more detail in appendix D. The notation is intended to emphasize the similarities 



^In general curved space-times questions of principle regarding the consistent definition of particle 
and composite-particle concepts ||3^ can arise. None of these affect our metric, however, due to its 
time-like Killing vector. 
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with the standard quantum numbers in flat space. The time coordinates of the center of 
mass are 

= ^(:^° + y°), ^ = 1,2. (46) 

The time dependence on the right-hand side of eq. (^) is entirely determined by time 
translation symmetry. For the same reason depends only on time differences, which 
we will, however, not exhibit explicitly. Finally, the normalization of the Bethe-Salpeter 
amplitude agrees with the normalization of the one-pion states adopted in appendix D. 

For the course of the following derivation it is convenient to introduce a finite pion 
mass, which allows the specialization to the pion rest frame. Later we will go back to the 
chiral limit. The bound state (pion) contribution to S*^^-*, obtained by inserting on-shell 
pion intermediate states into eq. (|4^) , can now be written in terms of the Bethe-Salpeter 
amplitude and its conjugate: 

. /• C?Po ^ X"(Xi,yi;p)»^x"(a^2,y2;p)^7 ,-^po(X?-X») /.yx 

After inserting this contribution into eq. (^) and specializing the external momentum to 
the pion mass shell, we end up with the homogeneous Bethe-Salpeter equation for x""- 

e~'^°^°x'^ixi,yi;p)af3 = J dfi{zi) J diJ,{z2) j dfi{z3) J dn^z^) 

In order to proceed further, the interaction kernel K has to be specified. The standard 
ladder approximation to the Bethe-Salpeter equation is consistent with the Hartree-Fock 
approximation for the quark propagator from the last section. We thus adopt it here and 
choose K accordingly to contain just the tree-level contact interaction of the NJL model: 

K{z,, Z2; zs, z,U,s^ = 2ig{Ti^T^, - Ti,T^^)6\z^, Z2)S\z2, Zs)6\z3, z,). (49) 

As a consequence, eq. (HSI) reduces (from now on we use matrix notation in fiavor, color and 



spinor space) to 
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-S{x, z)TWz, z- p)V^S{z, y)}e-'P'^\ (50) 
In flat space this equation has been solved analytically by Nambu and Jona-Lasinio in 



their original paper ||29|. The analytical solution ows its existence to the zero range of the 
interaction, and an analogous solution can thus be expected in the curved unit cell. 

We will now derive this solution explicitly from an infinitesimal chiral transformation of 
the quark propagator. Writing the chiral generators as = and using the explicit 

form of the self-energy (^) with the self-consistent Hartree-Fock propagator, we find 

Ts'^} = -2tg [ Ta tr [r^{5(^, z),T,^}] - {Siz, z),T^} Ta ] . (51) 

This transformation behavior follows directly from the chiral invariance of the NJL inter- 
action. (Note that the two coincidence limits t ±00 of S differ only in the vector part 
(cf. eq. (|26|) ). Since the latter drops out of eq. (^Tf), their further distinction becomes 
unnecessary.) Now we can write 



{S{x,y),T^} = J dfi{u) J df,{v)S{x,u){S-\u,v),T^}S{v,y) 
= -J dfi{z)S{x,z){^{z),T^}S{z,y) 

= 2ig j dii{z)S{x,z)[VAtA^A{S{z,z),T^}]~VA{S{z,z),T^}VA]S{z,y), (52) 

which has exactly the form of the Bethe-Salpeter equation ( [5D| ) for zero-momentum pions 
in the chiral limit, with the solution 

X"(x, y)=N {S{x, y),T^} = 2 AT j rf/i(z)5(x, z) S{z, y) (53) 

for the pion wave function. The normalization M will be fixed in the next section from 
the chiral Ward-Takahashi identity. Equation (^) remains a solution of the Bethe-Salpeter 
equation for a small, finite pion mass m,r, up to corrections of order m^. 

As the longest-wavelength excitations of the vacuum, the pions and their Bethe-Salpeter 
amplitude can be expected to depend strongly on the vacuum structure and its variations 
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with density. In the next section we will study this issue further by calculating the pion 
decay constant and its density dependence, building on the results derived above. 

VI. THE PION DECAY CONSTANT 

A second key parameter in the discussion of spontaneous chiral symmetry breaking, 
besides the quark condensate, is the decay constant of the pion, The explicit expression 
for the Bethe-Salpeter amplitude from the last section allows us to adapt the standard 
method of Nambu and Jona-Lasinio to calculate f-j^ in the curved unit cell. We start from 
its definition in terms of the axial current matrix element between the one-pion state and 
the vacuum, 

< 0|j5^(^)k'(p) >= -fJ^'d.riA^). (54) 

The space-time dependence of eq. (|5^) follows entirely from the spatial 5*0(4) symmetry of 
the hypersphere and from time translational invariance. It is contained in the eigenmodes r) 
of the Klein-Gordon equation on x R, which generalize the usual plane waves of Minkowski 
space and are explicitly given and normalized in appendix D, where also the set of generalized 
momenta p is defined. 

By covariantly differentiating eq. (Q) one obtains the hypersphere equivalent of the 
PCAC relation between the vacuum and one-pion states, 

< 0|V^j5",^(x)|7r^(p) > = -U6'''g--^d^g^^g,,d''v^{x) = Uml5^%{x). (55) 

In the derivation of eq. (|55D we used the free Klein-Gordon equation (|D2|) and the standard 
definition of the covariant derivative of a vector field, 

VM,(a:) = d>'A,{x) + r/''A^(x), (56) 

which contains the Christoffel connection introduced in eq. (^). The relation to the Laplace- 
Beltrami operator is immediately established with the help of the identity TJ^^^ = (d^g^). 
We further introduce the axial- vector and pseudoscalar current densities of the NJL model, 
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JsA^) = Qi^hf^l^Y^ix), j^ix) = g(x)i75yg(x), (57) 

which are connected by the divergence of the axial current, 

V''jl^{x) = 2mof,{x). (58) 

The relation between the pion decay constant and the Bethe-Salpeter amplitude, eq. (^), is 
established through the axial current matrix element (Q), which can be expressed in terms 
of the coincidence limit of x""- 

< 0\jUx)\Ap) >= -^==tT[^,l,-x''{x,x;p)]. (59) 

(The trace is over Dirac, flavor and color indices.) 

Before calculating from this relation, we have to fix the normalization of the pion wave 
function. To this end we consider the quark- ant iquark three-point functions of the current 
densities (|57D, 

Gl^{z,x,y) =< 0\T jl^{z)q{x)qiy)\0 >, (60) 

Gt{z,x,y) =< 0\Tr,iz)qix)q{y)\0 >, (61) 
and derive with the help of eq. (EBf) the chiral Ward-Takahashi identity 



Vf,) Gl^{z, X, y) = 2mo Gt{z, x, y) - S{z, y) 5\x, z) - S{x, z) 5\z, y) T^. (62) 

Integrating over z, the surface term on the left-hand side vanishes (recall that the pions have 
a finite mass) and we obtain 

2mo j dfi{z)Gt{z,x,y) = {S{x,y),T,-}. (63) 

The left-hand side receives contributions from one-pion intermediate states, which can be 
expressed in terms of the Bethe-Salpeter amplitude as 

2mo / dfx{z) Gl{z, x, y) = -iUx^i^, y), (64) 



28 



and, comparing eq. (^) with eqs. (|6^) and (|6^, we read off the normahzation constant 



(65) 



With the normahzation fixed, we can now equate the time components of the two ex- 



pressions for the axial current matrix element, eqs. (|5^) and (|59D , in the pion rest frame 
p = {0, 0, 0} with (cf. appendix D) 



-irtiTrt 



(66) 



to obtain an equation for the pion decay constant in terms of the Bethe-Salpeter amphtude: 



im-^fiT = -tr 

After inserting the explicit solution in the pion rest frame, eq. 



(67) 



we have 



-iniTrt 



tr 



7o75y5(0,x)75r3^(x,0) 



(68) 



We now use the Hartree-Fock Dirac propagator on 5^ x i? from section III to evaluate the 
right hand side further. The color and flavor traces in eq. ( p8| ) produce a trivial factor nji f 
and the complete trace (over color, flavor and Dirac indices) becomes 



tr[7o75y5(0,a;)75r35(x,0)] = 2n,nf [So{fi, -t)^2(/i,t) - So{fi,t)S2{fi, -t)] 

2i d 

= -4:ncnfSo{n,t)S2{iJ^,t) = —Tr:Siifi,t) 

m ot 



(69) 



(m is the dynamical quark mass), where we made use of the time-reversal properties of the 
propagator functions. 



S'o(/i, -t) = -Soifl, t), S'2(/i, -t) = S2{fl, t), 



(70) 



and of the identity 



— i 

S'o(/i,t) = — — S'2(/X,t), 

m ot 



(71) 



which follows directly from the explicit expressions for 5*0 and 5*2 in appendix B. After 
performing the two trivial angular integrations in the measure eq. (^^, and a partial 
integration in the time coordinate, the equation for the decay constant simplifies to 



29 



dte-'""-' d/i sinV ^2(/i,i)- (72) 

-oo JO 

We now use the spectral representation of the propagator function 5*2, eq. to write 

/oo r2n 
dte-'""^' dfx sin^ fx S^{ix,t) (73) 
-oo JO 

= 777 H / ^/^ *'^(/^) *™(^) / ' C^^) 

V4\// n^^iJO J-oo UJnUJm 



where 

Sn(Ai) = 2A;„L sm{knLfi) — cos(A;„L/i) tan ^. (75) 
The integral over /i can be done explicitly, 

/■27r 

/ c?/i s„(/i) s^(/i) = 27r(5„„L'„, (76) 
Jo 



with the quark level degeneracies D„ given in eq. (pQl). It remains to combine the above 
equations, to regularize the remaining mode sum consistently with the gap equation (by using 
the same regulator, eq. (p7D) and to introduce the filling factors of section IV to account for 
the valence quarks. Our final expression for the pion decay constant then becomes 

^^=|^i: ^"^';^"^ e.(A-fc„). (77) 

^ ^ n=l 

A comparison of this result with the fiat-space expression in ref. shows similarities 
analogous to those between the expressions for the quark condensate. In particular, 
gets contributions from the quark Dirac sea with the same functional dependence as in [Q. 
The main difference to the fiat-space result is again due to the new quark spectrum, which 
replaces the plane wave spectrum of euclidean space. 

We will continue the discussion of eq. ( |77D and its density dependence in section VIII. 
There we will also present results from the numerical evaluation and a comparison with the 
corresponding results of the chemical-potential approach. 



30 



VII. THE ZERO-DENSITY LIMIT 



At low baryon densities the volume of the unit cell becomes large and its curvature goes 
to zero. In the low-density limit, or equivalently in the large-L limit, one therefore expects 
to recover the results of the flat-space NJL model. In the present section we will perform 
this limit explicitly, both as a consistency check on our results and to illustrate how the 
quark spectrum and its degeneracies turn into the plane wave form. 

Let us start by considering the L ^ oo limit of the scalar self-energy, eq. (^3]). In this 
limit the Fermi momentum goes to zeroQ and contributions from quark propagation around 
the whole sphere are suppressed. Therefore it is convenient to use the representation of the 
propagator in terms of a sum over paths. Up to contributions from complete circles around 
the hypersphere we obtain 

— nm 1 — Ak'^T? 
= t4t^(1 + 2^^^/) / ©^(^ - ^) (78) 

^ ^(1 + 2n^nf) r dk — e,(A - k), (79) 

TT^ JO UJ 



which is indeed the original flat space result |29[. Similarly, the vector part of the self-energy 
becomes 

ZTT 1j Jo 

^ 0, (81) 

and goes, as expected, to zero in the large-L limit. 

To check the zero-density limit of is a bit more subtle, since the corresponding loop 
contains two quark propagators instead of the coincidence limit of only one. We start from 
eq. d?^ ) and use the representation ( |B23D for the propagator function 82- Neglecting again 
additional complete turns around the hypersphere, we have 



Recall that the Fermi momentum sets the baryon number in the cell to one. 
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POO r2n n 

fl = -Amncnf^amfL^ I dt ci/i[2/g(/i, t) + tan ^/o(/i, t)] 

Jo Jo 2 



incTifm 



2 " dko dk dq /•2^ 



47rL 



27r i-A 2tx j-a 277 io ^ 





2ikL + tan f 




2igL + tan | 


[kl - fc2 


- w? + ie] [/cq 




g2 _ _|_ ^gj 



(82) 



In the second line the time integration, which leads to a delta function 5{ko — go), and the 
consecutive go integration are already performed. Note that the range of the integration 
in the above expression can be changed to [— 27r, 0] without changing the integral, if we 
simultaneously change the signs of the integration variables k and g. We can thus replace 
half of the //-integral above by an integral with the alternative, equivalent range. After 
furthermore scaling // to /i = /iL, the L-dependent part of the above expression becomes 



^r^dkr^dq 

2L J-A 2tX J-A 27r J-2nL 

which, in the large-L limit, reduces to 

'•A dk dqiT 



2ikL + tan 



/i 



2iqL + tan 



/i 



(83) 



(84) 



Inserting this in the expression for above, we finally get 



J w 



—iricnfm 

ricnfm 
47r2 



A dk _ 
-A 2^ 

2 .A p 



dkn 



27r [/cq — /c^ — + ie] [kg — k"^ — w? + ie\ 



dk- 







(85) 



[uj = yW^+m?), which agrees with the flat-space result . 

The large-L limit thus reproduces, as expected, the standard zero-density quantities of 
the NJL model. The transition of the physics from the high- to the low-density regime can 
be traced to the differences in the quark propagation in small and large hyperspheres. At 
high densities geodesic paths with any number of turns around the sphere contribute and 
their sum builds up to the equivalent sums over the quark spectrum on S*^ x (cf. appendix 
C). In the large-L limit, on the other hand, only the direct path is relevant and yields the 
standard flat-space propagator. 
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VIII. RESULTS 



Before embarking on the quantitative discussion of our results, we fix the couphng g and 
the cutoff A of the NJL model at g = 4.08 GeV~^ and A = 700 MeV. These values have been 



found in ref. [|2| to best reproduce the phenomenological values of the quark condensate 
and of the pion decay constant in the vacuum. We adopt this choice to allow for a direct 
comparison between our results and those of ref. 0], and to ensure a phenomenologically 
acceptable zero-density limit of both the quark condensate and of /vr. Nevertheless, given 
the simplicity of the minimal NJL model and the sensitivity to the values of A and g already 
noted in ref. [Q, our numerical results should be considered as quahtative. 

The above model parameters could in principle be density dependent. Since we are 
mainly interested in new, qualitative features of the hj^ersphere formulation, however, we 
will follow the standard treatment and keep them constant. Finally, the diffuseness param- 
eter in the regulator, eq. (]27|), is set to e = 2.4. This choice compromises between a rapid 
convergence of the mode sums and a sufficiently smooth behavior of the observables at high 
densities. 

We are now ready to discuss our numerical results. The non-trivial solution of the 
gap equation (p7D, i.e. the dynamical quark mass m, and the negative cube root of the 
corresponding quark condensate, — < qq >^^^, are shown as a function of density in fig. 1. 
Both quantities decrease monotonically with density and finally go to zero at the density pc- 

The vanishing quark condensate does of course not necessarily imply a second-orderpha.se 
transition to chiral restoration. Recall that the condensate shown in fig.l was calculated by 
using the nontrivial solution of the gap equation in the quark propagator. It is not a priori 
clear, however, that this solution remains physical, i.e. the one of lowest free energy, if the 
baryon density increases. 



Indeed, it was claimed in ref. that in a certain parameter range of the standard 
NJL model, close to the chiral limit, the trivial solution becomes the global minimum of the 
free energy before the nontrivial solution reaches zero. This result, obtained in the chemical- 
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potential approach, would imply a first order phase transition to chiral restoration. However, 
studies in more realistic versions of the NJL model found a second-order transition for the 
whole phenomenologically acceptable parameter range |f44|| . 

To determine the order of the phase transition in our approach, we monitor the extrema 
of the free energy density 

^^^^ = ^(1 + 2n,n,) < qq < q^q £ - /.)u;„e,(A - K) (86) 

as a function of the cell volume or, equivalently, of the baryon density. The last term in 
eq. (|86|) is just the regulated sum over the energies of the occupied quark states, while 
the first two terms are required to avoid double-counting of the vacuum energy |^5|. It is 
straightforward to check that the extrema of eq. (^) under variations of the quark mass 
are , as expected, the solutions of the gap equation (|37|). 

The dependence of fi/V on the quark mass is plotted in fig. 2 for four baryon densities 
between p = and p = pc- These curves demonstrate that only two solutions of the gap 
equation exist at all densities up to pc- the trivial one at m = and the nontrivial one with 
finite m. In particular, the nontrivial (trivial) solution remains the minimum (maximum) 
of the free energy density up to Pc, where both extrema merge. The phase transition is 
therefore of second order. 

Closer inspection of the free energy density reveals that this result is, in contrast to the 
chemical potential approach, independent of the model parameters. Indeed, it is straight- 
forward to show that, as long as the nontrivial solution exists, it remains the minimum of 
the free energy density. This robust second-order phase transition, which actually persists 
for a much larger class of NJL models, is a welcome result, since the specific choice of the 
lagrangian is somewhat arbitrary. It also insures that a potential density dependence of the 
model parameters would not affect the order of the phase transition. 

For comparison, the quark condensate as given in ref. |^, i.e. calculated on the basis 
of the nontrivial solution in the chemical-potential approach, is also shown in fig. 1. There 
may exist parameter regions in which the transition of ref. [Q becomes actually first order 
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in the chiral limit. Introducing a current quark mass would in this case restore the second 
order transition without significant effects on the curve and the approximate critical density. 
In order to compare the density dependence and the critical densities of the second-order 
transitions, we therefore reproduce the original plots of ref. in fig. 1. 

Both approaches yield a very similar behavior of the order parameter with density. This is 
reminiscent of the analogous situation in the Skyrme model, where close similarities between 
hypersphere and array results first demonstrated the use of hypersphere calculations. 

In fig. 3 we plot the density dependence of the pion decay constant and again, for 
comparison, the same quantity as obtained in the chemical potential approach, f.^ decreases 
with density and goes to zero at the chiral restoration density. This is expected, since the 
pion has to decouple when the chiral condensate disappears. In the hyperspherical cell 
decreases initially somewhat faster then in the chemical potential approach, but goes to zero 
at almost the same critical density. 

Figure 4 shows the density dependence of both the quark condensate and the pion decay 
constant in the absence of valence quarks in the cell. This situation corresponds to setting 
/ci? = or, equivalently, to setting all /„ = 0. Even in this case both the condensate and the 
decay constant decrease and eventually vanish. This happens, however, at a considerably 
higher density. The absence of baryon number sources inside the cell delays the transition 
to chiral restoration, as one would expect. 

The valence-quark free cell can be interpreted as a region of reduced baryon density in 
the nuclear medium. While the interaction with the surrounding matter is still mediated by 
the shape and curvature of the 5 = cell, it is felt only by the quark vacuum in the interior. 
Such 5 = cavities can therefore neither be studied in the chemical-potential approach to 
the NJL model nor in the Skyrme model. In the latter, a unit of baryon number in every cell 
is crucial for the chirally restored phase to exist, at least in semiclassical approximation. In 
the standard approach to the NJL model at finite density, on the other hand, a homogeneous 
baryon density distribution is built in from the beginning. 

Due to the fixed relation between the ambient baryon density and the radius of all cells, 
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however, the size of the B — cavity cannot be varied at a given density. Such baryon- 
density free bubbles with delayed chiral restoration will therefore very likely not describe an 
equilibrium situation, but they might be related to a transient bubble phase in heavy-ion 
reactions, where the chirally broken and restored phases coexist. 
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IX. SUMMARY AND CONCLUSIONS 



In this paper we have generahzed a new approach to dense nuclear matter in the frame- 
work of extended hadron models. Originally developed for the Skyrme model, it divides the 
nuclear medium into cells of unit baryon number, which have the form of a 3-dimensional 
hypersphere S^. The new and crucial feature of these cells is their intrinsic curvature, which 
mediates interactions with the ambient matter. 

In order to explore the use of this approach beyond the Skyrme model, we first extend it 
to quark-based baryon models in general and then apply it specifically to the Nambu-Jona- 
Lasinio model, by investigating vacuum, constituent-quark and pion properties as well as 
their density dependence. The presence of fermions in the curved cell requires the introduc- 
tion of additional concepts, under which a new, isoscalar gauge interaction for the quarks is 
the most important. 

We find a high-density phase transition, which restores the chiral symmetry of the vac- 
uum by decondensing quark- antiquark pairs. At the same critical density the pions decouple, 
as signaled by the vanishing of the pion decay constant, which monotonically decreases with 
density. The phase transition is, for a large class of generalized NJL models, robustly of 
second order and independent of the model parameters. Furthermore, delayed chiral restora- 
tion takes place even inside of valence-quark-free cavities in the dense medium. In this case 
it is solely driven by the interaction of the quark vacuum with the cell curvature. 

Remarkably, we find close similarities between the density dependence of our results and 
those of the conventional chemical-potential approach. The behavior of the constituent quark 
mass and of the quark condensate with density, for example, is (in the phenomenologically 
reasonable NJL parameter range) in both cases almost identical. This similarity may seem 
at first surprising, since the two approaches are based on rather different descriptions of the 
dense environment. At least qualitatively, however, it is easy to understand. Recall that the 
energy gain from the Dirac sea, filled with massive as compared to massless quarks, leads 
at zero baryon density and for sufficiently large coupling to the development of a nontrivial 
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vacuum. In the chemical potential approach, this energy gain gets with increasing density 
more and more compensated by the occupied states in the Fermi sea, since they count with 
opposite sign in the spectral sum of the gap equation. At the critical density, it becomes 
energetically favorable for the quarks to be massless, and the chiral phase transition occurs. 

On the hypersphere the mechanism is different in principle, but similar in effect. The 
occupation numbers of the valence quark levels are here density-independent, but the spec- 
trum itself changes in a characteristic way. Due to the stronger localization in the smaller 
high-density cells, the quark momenta grow inversely with the cell size. This leads to a 
reduced number of states in the regularized Dirac sea, with the same consequences as in the 
chemical potential approach. 

The overall increase of the quark momenta at larger densities is an effect common to all 
compact unit cells. The specific and detailed change of the quark spectrum, however, and 
therefore the quantitative density dependence of the results, is governed by the particular 
form of the unit cell. The close agreement with the chemical-potential approach thus sup- 
ports the choice of as the cell geometry. An analogous agreement between results of the 
conventional and the hypersphere approach was observed in the Skyrme model. 

Experience gained from hypersphere calculations in different models can help to disen- 
tangle features specific to the Skyrme model from more general or even model-independent 
aspects of the hypersphere approach. In this respect our results indicate, for example, that 
the presence of a soliton is not required for the hypersphere approach to be effective. In 
particular, the specific topological properties of the Skyrme soliton and its hedgehog form 
are not indispensable. Since the winding number of the skyrmion plays an important role 
in its hypersphere behavior, this was not a priori obvious. Furthermore, a non-zero baryon 
number [i.e. winding number) in the cell and the corresponding hedgehog structure are 
necessary requirements for chiral restoration in the Skyrme model, whereas this condition 
is suspended in the NJL model, as the results in the B — cell show. The NJL results 
also indicate that the hypersphere approach is compatible with different mechanisms for 
spontaneous chiral symmetry breaking. 
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The application of the hypersphere approach has often practical advantages, which may, 
however, vary in different models. In the Skyrme model, for example, 5*^ calculations re- 
quire drastically less numerical effort than conventional unit-cell calculations. While this 
particular benefit does not translate to the NJL case (since the numerical requirements of 
the conventional (chemical-potential) approach are already moderate), the hypersphere ap- 
proach to the NJL model can be applied in realms not easily accessible by other methods. 
The discussed low-density bubbles in nuclear matter with delayed chiral restoration can e.g. 
neither be studied in the Skyrme model nor in the chemical-potential approach to the NJL 
model. Another profitable application of the hypersphere approach to the NJL model is 
encountered in the study of nuclear matter at finite temperature. While the Skyrme model 
calculation meets with rather involved technical difficulties in this situation even on the 
hypersphere, it can be straightforwardly addressed in the NJL case . 



We regard the above results as promising indications for the hypersphere approach to 
be useful beyond the Skyrme model. Nevertheless, many important and interesting ques- 
tions, also on a more fundamental level, remain open. It would, for example, be useful to 
know if hyperspherical and flat unit cells can be explicitly related to each other, at least 
approximately. A hint in this direction comes from the work of ref. [^, where matter at a 
fixed baryon density is almost identically described in hyperspherical cells with both baryon 
number and volume doubled, as a first step towards the limit of a skyrmion matter configu- 
ration in fiat space. One might also hope to make contact with older speculations ||52| about 
a description of chiral dynamics in terms of space curvature. 

The theoretical framework developed in the present paper can straightforwardly be 
adapted to other chirally symmetric quark and hybrid models, as for example to bag models 
|53| and to non-topological soliton models [Q. It might also be interesting to study the 



relations between the Skyrme and NJL models, which can be established by bosonization 



of the latter |3^,^^, in the hypersphere formulation. Finally, our work gives an explicit 
example for symmetry breaking and restoration due to space-time curvature in a fermionic 
theory, and as such may be also useful in other areas of physics. Related questions are, for 
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example, presently discussed in a cosmological context []55 
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APPENDIX A: MAURER-CARTAN BASIS AND DIRAC SPECTRUM ON x R 



The generalization of Lorentz (Poincare) tensor fields to curved space-times is straight- 
forward. One simply replaces the Lorentz group by the group of general coordinate trans- 
formations, which becomes locally {i.e. in the tangent spaces of the manifold) the full linear 
group GL{R, 4). The generalized tensor fields form its representations. 

The situation is more involved for fields of half-integer spin, because the linear group 
has no double- valued {i.e. spinor) representations. For this reason, spinors have to be 
introduced with respect to a local, orthonormal basis, the vierbein (in the (co)tangent 
space of space-time). As discussed in the introduction, the metric and therefore the physics 
is independent of the arbitrary orientation of these frames, which can be locally changed 
by a Lorentz gauge group. The fermion fields form the spinor representations of this local 
Lorentz group, which acts in its defining representation on the Latin index of the vielbeinQ. 

For explicit calculations we have to choose an appropriate vielbein on x R, i.e. we 
have to fix the gauge of the local Lorentz group. The spatial hypersphere has the geometry of 
the group manifold of SU{2), scaled to radius L. This suggests the rescaled Maurer-Cartan 
vector fields of SU(2) as the natural choice for the spatial components of the vielbeinQ 

ea = ^TT[Ta{d,h)h-'] g'^ d, {a,t = 1,2,3), (Al) 

where we parametrized the elements of SU{2) in polar coordinates as 

h{x) = cos /i + iTaf"'{0, 0) sin fj,. (A2) 



This local definition of the spinors can in general not be consistently extended to the manifold 
as a whole. x R, however, as a group manifold, is parallelizable. In technical terms this implies 
the vanishing of its first two Stiefel- Whitney classes, which guarantees the existence of a global spin 
structure. 

^^To be specific, we use the right-invariant Maurer-Cartan fields. This is purely a matter of 
convention. 
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{f°'{6, 0) is the unit vector in R^.) Complemented by the time component Co = dt we obtain 
an orthonormal vierbein on S'^{L) x R, from which the metric eq. (|1|) can be recovered 
with the help of eq. (^. The simplest way to derive the corresponding Levi-Civita spin 
connection is to solve Cartan's first structure equation ||3^ in the same "Maurer-Cartan 
gauge"Q. The result is 

uJa = j^GabcCr'"' (o, 6, c e {1, 2, 3}) and Uo = 0. (A3) 

{o'ab = f[7a,76]) Notc that the Ua are space-time independent and essentially given by the 
structure constants of SU{2), which is the main advantage of using the Maurer-Cartan 
vielbein. 

The Dirac operator on x R, which appears in the free part of the NJL lagrangian, eq. 
(p!^), and to which we add a constant vector self-energy (we do not include the scalar part 
of the self-energy, since it will later be absorbed in the mass term), is 

TJp=[ir{ea + UJa) + Klo] (A4) 

and now takes the simple formQ 



lip 



I A.^ 



A_ 



(A5) 



^^A fixed vierbein does in general not determine the corresponding spin connection completely. We 
take our cell to be torsion-free, however, and further require the spin connection to be compatible 
with the metric (metricity condition), so that it is uniquely given by the first Cartan structure 
equation. 

^^We use the chiral representation of the 7-matrices: 7° = /? = , 7 = , a 
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in terms of the operators 



(A6) 



(Arrows indicate vectors in the spatial, euchdean subspace of the local vielbein.) 

For the interpretation of our results it will be useful to derive the energy spectrum of 
the corresponding Dirac equation. Let us therefore look for the static solutions 



-iuint 



(AT) 



of the Dirac equation 

/ 

—m A_|_ 
A_ —m 

which has the equivalent hamiltonian form 



{Hp - m)^r, 



\ 



( ^(1) \ 



0, 



(A8) 



[-ia{e - uj) + (5m\ i/jni^) =a;„'0n(x), 



(A9) 



where a;„ = cDn + E^,. (The constant vector self-energy just shifts the origin of the energy 
scale.) Now we can again exploit the properties of S^{1) as the group manifold of SU{2). 
The invariant vector fields of S^{L) are (as for every Lie group) simply the group generators 
La, rescaled to radius L: 

2 



Ca = —La 

iL 



(AlO) 



where the !/„ satisfy the usual angular momentum commutation relations 



[Lai Lh\ — i^abcLc- 



(All) 



The appearance of a "spin-orbit" term in 



A^r = {-'^ -K± ]r{4SL + h}r a = 1,2 



(A12) 



[S is the spin operator (t/2) then suggests to diagonalize A-t by taking in the coupled 
basis with J = L-\- S, so that 
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A' 

7=/ + i. / = 10.1.2...! 

(A13) 



\l r,,^ for j=/ + |, / = {0,1,2...} 



+ V^l,^ for j = /-i, / = {1,2,3...}. 
Inserting eq. ( [A13|) into the Dirac equation and combining the resulting two coupled equa- 



tions for the upper and lower spinor components one obtains 



(2/ + 1)2 for 3 = l + \ 



< = m' + T^^{ ' ' (A14) 



L2 



-21 - \f for ]=l-\ 
which determines the energy eigenvalues uOn- All properties of the spectrum can be read 
off directly from eq. ( |A14 ): each pair of levels with j = ^,1 = and j = = f. 



n = {1,2,...} contains 

Dn = 2n{n + 1) (A15) 

degenerate states with momentum fc„ and energy iVn, 

kn = ^l = kl + m\ (A16) 

Finally, the total number of states up to the A^th level sums up to 

N 2 

Y,Dn = -N{N^ + ?,N + 2). (A17) 
1 

This spectrum agrees with existing derivations in the literature, obtained by various different 



methods, see e.g. [^] and also |^ for the massless case. Let us finally mention that the 
existence of a time-like Killing vector (cq = dt) in our cell permits a unique definition of in- 
and out-states (and their Fock spaces) on the basis of this spectrum. 

APPENDIX B: THE FERMION PROPAGATOR x R 

In this appendix we construct the quark propagator S{x) on 5*^ x i? for a Dirac operator 
of the form discussed in appendix A, i.e. including a constant self-energy ([5^). S{x) is the 
solution of 
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{Ul)-m)S{x,y) = 5\x,y), (Bl) 
with Up given in eq. (|A4|). The scalar part of E is combined with the current quark mass 



in the mass term m = ttlq + S^, the vector part is absorbed in the covariant derivative 
and Feynman boundary-conditions are imphed. (We will later generalize the boundary 
conditions to take the presence of valence quarks into account.) 

The four-dimensional delta function is generalized to curved space by demanding 

Jdf,{x)6\x,y)fix) = f{y), (B2) 

so that 

S\x,y)^{-g)-'/'l[S{x'^-y>^), (B3) 

11=0 

up to the sum over periodic paths (see below), (g = — sin^ /x sin^ 6* is the determinant of 
the metric.) 

Due to the compactness of the unit cell, the free propagation of a quark between x and y 
can proceed either directly (on a geodesic) or via an arbitrary number of intermediate turns 
around the hypersphere. We can therefore consider the geodesic coordinates as non-periodic 
real numbers and count the number of windings of a given path in multiples of 2tt. The 
propagator then becomes a sum of the propagators for each of these paths. (In order not 
to complicate the following equations unnecessarily, we will suppress these periodic sums 
until they play an active role in the derivation.) The above procedure ensures that the full 
propagator is periodic in the coordinates, as required by the physical situation, and allows to 
write its Fourier transform in terms of continuous momenta. Although the Fourier integrals 
can be transformed into equivalent sums over the discrete quark spectrum in the compact 
space (see appendix C), the former representation turns out to be more convenient for the 
derivation of the propagator. 

To derive the Dirac propagator explicitly, we start from the representation 

S{x,y) = -{^ + m)G{x,y), (B4) 
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which defines the block-diagonal 4x4 matrix 

/ 

G{x,y) = 



V 



(B5) 



Gix,y) 
G{x,y) 

a Greens function of the iterated Dirac operator. By inserting eq. ( [B4| ) into the defining 
equation (|B1| ) for S{x,y), we then obtain 



-A+A_ + m')G{x,y)=S\x,y), 



(B6) 



with A-i- given in eq. ( [A 61 ), and 



A+A_ = (idt + + ee- -ae 



4L2' 



(B7) 



(The derivatives act on x.) Note that ee is the Laplace-Beltrami operator on 5*^: 



ee = g 2 dig2 g'^ dj 

92 d 1 

— h 2 cot fi— — I — 



(B8) 



1 



92 



dfi ' sin2 fj, \ 86'^ 86 I sin^ /i sin2 9 dcjP ' 



The 2x2 matrix structure of G{x, y) can be expanded in the Pauli-matrix basis, 



G(x, y) = Gi{x, y) + aeG2{x, y), 



(B9) 



where Gi,G2 are scalar functions. 

As in Minkowski space, the required invariance of the vacuum under the symmetries of 
the metric leads to considerable simplifications in the space-time dependence of the prop- 
agator. Due to the symmetry group of the spatial hypersphere, 5*0(4), free propagation 
between two points is equivalent to propagation of the same geodesic distance from an ar- 
bitrarily fixed pole0. Taking the pole as the origin of the polar coordinate system on S^, 



^^This can be made explicit with the help of the seven Killing vectors of x R, which generate 
the isometry group. They replace the ten Killing vectors of the Poincare group in flat Minkowski 
space. 
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the geodesic distance simply becomes /iL. Furthermore, the propagator will depend only on 
time differences t (the reference time is chosen to be zero), due to the invariance of both 
metric and dynamics under time translations. 

After inserting the decomposition ( P9| ) of G into eq. ([B6D, we obtain two coupled equa- 
tions for Gi^2- 

9 



and 



-{idt + - e-e + — + m^) t) + '-ee G2{^^, t) = 6\x) (BIO) 



{idt + -^+^ + ^2(/U, t) + ^ Giifi, t) = 0. (Bll) 



(In the above derivation we used the Lie bracket relation [ea,efe] = j^e^j^ec, see eqs. ( |A10D , 



( |A11|) .) The remaining 6 and (p dependence in eq. (pi Op is removed by integrating the first 
one over both angles with their corresponding measure. (Eq. ( [Bll| ) is already independent 
of 6 and 0). The remaining explicit /i dependence can then be absorbed into the functions 
G'i,2(/^5 1) = sin yU Gi^2{f^, i) by writing the relevant part of ee as sin""*^ yu(H-^) sin fi and 
the /i-dependent part of the delta function as — sin^^ fi6' (fi) (a prime indicates a derivative 
with respect to fi). The Fourier transforms of 6*1.2 then satisfy simple algebraic equations, 

(fco + E,)2 -e-m'--^) G,{ko, k) + y (e - 4) G2{ko, k) 



((fco + S.)' - - + iZ^) ^2(^:0, k)-'- Gi{ko, k) = 0, (B12) 
which have the solutions (we introduce the abbreviations k± = k± 6 = {ko + S^)^ — m^) 



and 

G,{ko, k) = -iL {5-k^ + ^) G2{k,, k). (B14) 
After transforming back to space-time and imposing Feynman boundary conditions, we get 
Gi(/i, ^) = ^ (2 cos |r(/i, t) + sin ^/(/z, t)) e^^"* (B15) 
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and 



G2(/i,t) = -^sin^e^^''*J(^,t) (B16) 

in terms of the integral^ 

which is the scalar Feynman propagator in two-dimensional Minkowski space. The ultra- 
violet regularization of this integral and the change in the pole structure of the integrand 
due to valence quarks will be considered in appendix C. The full (2 x 2) propagator can now 
be obtained by inserting eqs. (P15|) and (P16|) (divided by sin //) into eq. ( pOf ) and carrying 
out the periodic sum: 

+00 

G{x) = a(sin/i)-i ^K^Sf+^.t) ^ {-iy[2I'^ + tan ^4], (B18) 

n=— 00 ^ 

where In{fJ',t) = + 2mT,t), the prime represents a derivative with respect to /i and 
a = I/AtcL"^. It remains to evaluate the expression ( P4|) for the spinor propagator, using the 
block anti-diagonal form of the Dirac operator given in eq. (|A5|) . As an intermediate step, 
we calculate 



-00 



A±G'(x) =a(sin/i)-^e'(t-^+S''*) ^ (-l)"[2j; + tan 



T^(2/:-cot|/0], (B19) 



where the dot represents a time derivative. Inserted into ( |B4| ) , we obtain the quark prop- 
agator explicitly. It can be brought into a familiar form by expanding the Dirac-matrix 
structure in the 7-matrix basis: 

S{x) = e^(^^^+^"*)[5o(/x,t)7o - Si{fi,t)n-S2ifi,t)]. (B20) 

(E is the Dirac spin matrix with the Pauli matrices on the diagonal and should be distin- 
guished from the self-energy.) The three invariant scalar amplitudes Si are defined as 



^^lifj-, t) can be expressed in terms of a Hankel function. 
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^o(/i,t) = -^«(sm^)-i 5: (-ir[2/; + tan|jJ, (B21) 

n=—oo 

S^M=^^{sin^,)-' £ (-1)"[2/: - cot ^J^], (B22) 



n=— oo 
+00 



S2ifi,t)= ma {sin fx)-' ^ (-l)"[2/; + tan|4]. (B23) 



Eqs. ( |B2U| ) - ( [B23D are the central result of this appendix. The closed, analytical form of 



the fermion propagator on x R may be useful also beyond the context of the present 
paper. 

The sums over all geodesic paths connecting the pole with x ensure, as mentioned above, 
the required spatial periodicity of the propagator. This can now be easily checked explicitly. 
The J„(/x, t) and their derivatives, appearing in eqs. ( P21|) - (|B23| ), satisfy 

J„(^ + 27r,t) = J„+i(/i,t). (B24) 

Inserted into the expression for the quark propagator, the expected periodicity in /i follows 
immediately: 

S{ii + 2n, 9, 0, t) = (-1) e''^^' S{fi, 9, 0, t) = 9, 0, t). (B25) 



APPENDIX C: SPECTRAL REPRESENTATION AND COINCIDENCE LIMITS 

The sums over geodesic paths in the propagator functions (P21|) - (|B23|) can be rewritten 
alternatively as sums over the spectrum of the Dirac operator. For many applications the 
spectral representation, which we will now derive, is more transparent and convenient. In 
numerical calculations it is particularly useful at high densities, where the cell is small and 
paths which circle many times around the sphere contribute significantly. 

The n dependence of the integrals /„, defined after eq. ( |B18|) , and their derivatives 
resides exclusively in the factors (—1)"" exp{2i7mkL). Exchanging the order of the sums over 
n and the integrals over k and using the identity 
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L E (-1) 



n^iirnkL 



E 5(k- 



2m +1 
2L 



E S{k-km), 



(CI) 



we can replace the sums over geodesic paths in the propagator by mode sums over the 
discrete quark momenta, eq. ( |A16| ). After performing the now trivial k integrals this leads 
to the spectral representation. Recall that the integrals over continuous k in the compact 
cell appeared since we took /i effectively in the interval fi G [0, oo], which allowed the use of 
continuous Fourier methods in the derivation of the propagator. Consequently, we had to 
add contributions from any number of (physically indistinguishable) geodesic circles. The 
two equivalent ways of representing the compactness of the unit cell are simply related by 
eq. dCl. 

For the propagator functions 5*0 and 5*2 we now get the spectral representations 



ma 
2^ 



+ 00 „ 

(sin/i)"^ E J 



dkn 



J{k„LiJ,-kot) 



2ti kl 



m^ + it 



2iknL + tan 
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2knLsm{knLn) — cos(/c„L/x) tan ^ 
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(C2) 



and 



S'o(/i,t) = 9t S'2(/i,t) 
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E(sin/x) ^ 



n=l 



2knL sin{knLfi) — cos(A;„Lyu) tan 



/i 



X 



0(t)e 



-iuJnt 
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(C3) 



The evaluation of the constituent quark self-energies in section IV requires the regularized 
coincidence limits of S. The time-ordering ambiguity of the propagator in the coincidence 
limit is resolved as usual by referring to the order in which the fields appear in the interaction 
lagrangian. Accordingly, we define 



S^{0)= lim ^(xo,x = 0) = 50^70-52 

a;o— »0± 



and obtain 



^0 = TTT E(Tl + /n)^ne,(A - kn), 
n=l 



s. 



(C5) 
(C6) 
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Here we have additionally accounted for the presence of the valence quarks by adding the 
filling factor terms introduced in section IV. We recognize again the appearance of the quark 
spectral properties, which can alternatively be derived directly from the Dirac equation (see 
appendix A). 

APPENDIX D: THE FREE SPIN-0 FIELD ON x R 

As a prerequisite for our discussion in sections V and VI, we will derive here the mode 
decomposition of the free pion field (or any real spin-0 field, in general) in the hyperspherical 
unit cell. We denote the set of three spatial coordinates by (x) and the three generalized 
momentum quantum numbers (see below) by (k), and we suppress the trivial isospin indices. 

Let us now expand the free field in the x R background as 

<Pix) = (^k(a;)ak + Vl{x)al) (Dl) 

k 

in terms of the complete set of solutions r^k of the free Klein-Gordon equation]^ on x R, 
which generalize the usual plane waves of Minkowski space: 

(□ + m^) r]^{x) = - ee + m^^ r]^{x) = 0. (D2) 

(e e is the Laplace-Beltrami operator defined in eq. (B8).) The r/'s are chosen to be or- 
thonormal in the generalized scalar product 

{Vk,Vk') = -i J^^ dfi{x) r]Ux)—r]l.,{x) = (5k,k'- (D3) 

Note that the integral above does not include the time direction and is restricted to the 
spatial hypersphere. One easily shows that this definition of the scalar product is time 
independent. In order to find the explicit form of the modes ?7k, we first separate the time 
dependence. 



i^Note that we consider the minimally coupled case. In general an additional term proportional 
to the Ricci scalar of the metric can be added to the Klein-Gordon equation. 
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r]k{x) = Mk(x)-^==, (D4) 

and then solve the remaining static eigenvalue equation of the Laplace-Beltrami operator: 

eetik(x) = -K^nk(x) = (m^ - ul) Mk(x), (D5) 

where we defined ul_ = + m'^ . Separating further the dependence on the spatial coordi- 
nates, the general solution of eq. (|D5| ) can be found to be 

Mk(x) = , ^ sinV Ci'+^Hcos/i) YiUe,(j)). (D6) 



Here the C^") are ultra-spherical Gegenbauer polynomials the Yi^ are spherical har- 
monics and the h!^^ are normalization constants defined by 

n\[n + ajl ^{a) 

The set k is now specified by the three quantum numbers k = {n, I, m} within the range 

n = 0,l,2,... 
/ = 0,1,2,... 

m = -1,-1 + 1,.. .,1-1, 1. (D8) 

The eigenvalue corresponding to the eigenmode is 

4 = <i = ^{n + l){n + l + 2). (D9) 

The static modes are orthogonal on 5*^ and normalized to 

/ (i/x(x) Mk(x)Mk,(x) = 5k,k', (DIO) 

so that the time-dependent modes ?7k satisfy the orthonormality relation eq. ( p3| ) . The 
standard equal-time commutation relations 

[$(x,t),n(x',t)]=^53(x,x') (Dll) 
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(n = dt^ is the canonical momentum conjugate to 6^ is the spatial part of the delta 
function defined in appendix B) then insure the norm of the one-particle states: 

[ok, a|^/] = 5k,k' 

k, Ok'] = k, al] = 0, (D12) 

so that 

< 0\(j)''{x)\TT\k) >= (D13) 
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FIG. 1. — < qq >^/^ as a function of density a) in the hyper spherical unit cell (solid line) and 
for camparison b) in the chemical potential approach of ref. [Q] (dashed line). Also shown is the 
constituent quark mass (dotted line). 



FIG. 2. The free energy density (in units of W^'^MeV^) as a function of the quark mass for 
four values of the baryon density. 



FIG. 3. /tt as a function of density a) in the hyperspherical unit cell (solid line) and for 
comparison b) in the chemical potential approach of ref. @ (dotted line). 



FIG. 4. - < qq >^/^ and (in MeV) function of density in the B = cell, which 
contains no valence quarks. For comparison, the same curves in the B = 1 cell from fig.s 1 and 3 
are also shown (dotted line). 
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